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from  0 to- 6 degrees. 
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These  Unique  data  help  eo  fill  an  Important  gap  in  our  experimental 
knowledge  of  near-surface,  low  speed  hydrofoil  performance.  ' This  speed  range 
has  only  recently  attracted  attention  because  of  the  large  sizes  of  projected 
hydrofoil  lift  systems.  The  present /experiments  indicate  that  there  are 
significant  changes  In  lift  and  drag  for  operation  at  shallow  and  moderate 
submergences  in  the  low  chord  Froude  number  range,  especially  in  the  vic-lwti!y 
-of  0.75  to  1.5. 

A prediction  of  the  residual  drag  coefficients  for  an  18-foot-chord 
hydrofoil  operating  at  constant  lift  shows  that  the  wavemaking  drag  Is 
significant  and  that  the  trend  of  tire  Increase  In  the'  residual  drag  at  Froude 
numbers  ranging  from  3.5  down  to  1.75  appears  to  be  roughly  parallel  to  the 
trend  of  the  induced  drag  alone.  However,  the  magnitudes  of  the  residual 
drag  coefficients  are  rather  large,  and  there  are  Indications  that  Interference 
wave  drag,^due  to  the  presence  of  a strut  and  a pod,  are. an  important  feature  of 
the  low  speed  drag  variation. 
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ABSTRACT 


Results  of  towing  tank  experiments  are  presented  for 
lift,  drag,  moment,  and  center  of  pressure  of  an  aspect  ratio  4, 
rectangular  planform,  5.21  percent  thick,  flat  plate  hydrofoil 
attached  to  a large  strut  and  pod  and  operated  beneath  the  free 
surface  of  water.  The  nominal  range  of  the  chord  Froude  number 
was  0.5  to  3.5  with  corresponding  chord  Reynolds  numbers  of 
0.74  X 10^  to  5.2  X 10^.  Tests  Included  seven  foil  submergence 
depths  ranging  from  0.25  to  3.5  chord  lengths,  and  foil  angles 
7f  attack  varying  from  0 to  8 degrees. 

These  unique  data  help  to  fill  an  Important  gap  in  our 
experimental  knowledge  of  near-stirface , low  speed  hydrofoil 
performance.  This  speed  range  has  only  recently  attracted 
attention  because  of  the  large  sizes  of  projected  hydrofoil 
lift  systems.  The  present  experiments  Indicate  that  there 
are  significant  changes  in  lift  and  dra,.:  for  operation  at 
shallow  and  moderate  submergences  in  t^  ow  chord  Froude 
number  range,  especially  in  the  vicinity  • ' 't.75  to  1.5. 

A prediction  of  the  residual  drag  co.i f Iclents  for  an 
18-foot-chord  hydrofoil  operating  at  constant  lift  shows  that 
the  wavemaking  drag  is  significant  and  that  the  trend  of  the 
increase  in  the  residual  drag  at  Froude  numbers  ranging  from 
3-5  down  to  1.75  appear  to  be  roughly  parallel  to  the  trend  of 
Che  Induced  drag  alone.  However,  the  magnitudes  of  the  residual 
drag  coefficients  are  rather  large,  and  there  are  indications 
that  interference  wave  drag  due  to  the  presence  of  a strut  and 
a pod  are  an  Important  feature  of  the  low  speed  drag  variation. 
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ADMINISTRATIVE  INFORMATION 


This  work  was  sponsored  In  part  by  the  Advanced  Naval  Vehicles 
Concepts  Evaluation  (ANVCE)  Project,  Task  Area  SS  H15002,  Work  Units 
1-1861-014-4A  and  1-1102-003-44,  and  by  the  Naval  Sea  Systems  Command, 
Task  Area  SF  43  421  202,  Work  Unit  1-1500-100-13. 

INTRODUCTION 

As  the  displacement  and  size  of  proposed  subcavltating 
hydrofoil-supported  ships  have  grown  larger  and  larger  in  recent 
years,  the  physical  sizes  of  the  foils  required  for  these  designs 
have  naturally  also  grown  considerably.  For  typical  speeds  near 
takeoff,  at  around  20  to  25  knots,  the  corresponding  chord  Froude 
number  range  of  interest  has  shifted  to  values  lower  than  are  available 
in  the  existing  data  base  of  hydrofoil  hydrodynamic  characteristics. 
Chord  Froude  numbers  approaching  F^  = 2 or  1.5  are  now  not  unreasonable 
values  to  consider  as  extremes  for  takeoff  performance  of  very  large 
hydrofoil  systems. 

The  distinctive  features  of  the  hydrodynamic  forced  on  a hydrofoil 

traveling  near  a free  surface  at  low  Froude  numbers  are  a rapid  rise 

in  the  curve  of  induced  plus  wavemaking  drag  coefficient  accompanied  by 

a dip  in  the  lift  coefficient.  These  properties  had  been  revealfid 

long  ago  by  calculations  using  linearized  potential  theory  as  presented 

2 

for  example  by  Nlshlyama,  summarized  in  Reference  l,by  Wu  and  by 
3 

Breslln.  For  Instance,  Figure  1 displays  the  predicted  behavior 


2 


for  lift  and  total  Invlacld  drag  coefficients,  as  determined  In 

It 

computations  by  Nlshiyama  for  a rectangular  planform  wing  with  a 
foil  submergence  ratio  f/c  " 0.59.  In  these  plots,  the  C is  the 
resultant  lift  coefficient  for  the  foil  operating  In  the  presence  of 
the  free  surface,  is  the  reference  lift  coefficient  for  the  same 

oo 

foil  in  an  unbounded  flow,  and  (C^  + Is  the  total  Induced  drag 

(including  biplane  effect)  plus  the  wavemaking  drag  coefficient. 

The  general  trends  of  these  predicted  features  have  been  partially 

verified  experimentally  down  to  » 1 at  f/c  *■  0.84  for  an  aspect  ratio  10 

3 

hydrofoil  as  pointed  out  by  Breslln  in  compa»‘lng  his  theory  with  the 
results  of  tests  on  rectangular  planform  hydrofoils  presented  by  Wadlin, 
et  al.^  It  may  be  noted  that  the  low  speed  end  of  the  test  results 
reported  in  Reference  5 are  somewhat  suspect  because  the  section  drag 
coefficients  (total  viscous  drag)  are  clearly  smaller  than  the  turbulent 
flat  plate  friction  line.  Guaranteeing  that  the  measured  friction  drag 
always  falls  on  the  same  transition  curve  between  laminar  and  predominately 
turbulent  flow  is  a notoriously  troublesome  experimental  problem  for 
experiments  conducted  at  Reynolds  numbers  in  the  range  5 x 10^  to 
6 X 10  ^without  turbulence  stimulation. 

Most  hydrofoil  experiments  are  designed  to  cover  relatively  high 
Froude  numbers.  For  example,  the  extensive  hydrofoil  test  program 
reported  on  by  Feldman^  produced  data  for  the  combination  of  a foil 
plus  strut  (s)  for  rectangular  planform  hydrofoils  having  six  different 
aspect  ratios,  at  six  depths  of  submergence,  with  chord  Froude  numbers 
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ranging  between  2.06  and  9.63.  More  recently,  Layne^  has  performed 
experiments  on  tapered  planform  hydrofoils  similar  to  the  forward 
foil  of  the  PCH  having  two  different  NACA  section  shapes.  The  chord 
Froude  number  range  covered  was  7.44  to  23.2. 

As  far  as  is  known,  there  exists  no  systematic  data  for  hydrofoil 
performance  at  speeds  spanning  the  low  Froude  number  range  where  the 
most  dramatic  speed-dependent  changes  occur.  The  experiments  reported 
on  here  have  been  planned  to  provide  an  initial  body  of  Information 
specifically  for  the  low  Froude  number  speed  regime. 

DESCRIPTION  OF  MODEL  AND  APPARATUS 

HYDROFOIL 

The  hydrofoil  model  built  for  these  experiments  was  made 
Intentionally  simple  for  ease  of  manufacture  and  low  cost.  It 
consisted  of  a rectangular  planform,  flat  plate  of  6061-T6  aluminum 
having  a thickness  of  1.25  inches  (0.03175  m)  with  a 0.625  inch 
(0.0159  m)  leading  edge  radius  and  a wedge-shaped  afterportion 
tapered  to  a sharp  trailing  edge  having  an  included  angle  of 
21.8  degrees.  With  a span  of  8 feet  (2.438  m) , the  plate  had  simple 
squared-off  tips.  The  2 foot  (0.6096  m)  chord  length  was  chosen  to 
provide  the  desired  Froude  number  range  using  the  useful  speed 
capability  of  the  towing  carriage,  and  it  allowed  reasonable  Reynolds 
numbers  (with  turbulence  stimulation)  at  the  lowest  speeds,  while  at 
the  same  time  keeping  the  span  at  a manageable  size.  It  also  limited 


] 
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the  magnitude  of  the  largest  expected  lift  forces  to  values  that 
could  be  handled  with  available  block  gauges.  The  foil  surface 

\ 

had  a rolling  mill  finish,  and  the  regions  near  the  machined  leading 
and  trailing  edges  were  smoothed  with  emery  cloth.  Table  1 summarizes 
the  main  features  of  the  flat  plate  hydrofoil  geometry.  I 

' In  order  to  Insure  fully  turbulent  boundary  layer  flow  throughout  j 

i 

i the  speed  range,  wire  trips  were  placed  at  the  nose  of  the  foil, 

I 

! extending  from  tip  to  tip  across  the  entire  span  on  both  the  top  ; 

and  bottom  sides.  It  was  felt  that  the  turbulence  stimulators  should 

be  located  somewhat  forward  of  the  most  likely  position  of  flow 

separation  which  for  an  isolated  circular  cylinder  Is  known  to  occur 

8 

I at  about  81  degrees  along  the  arc  of  the  surface.  The  two  trip 


wires  were  placed  nominally  at  + 60  degrees  along  the  nose  radius, 
or  at  a distance  of  0.312  inch  (0.00794  m)  back  from  the  leading  edge. 
The  standard  criterion  for  sizing  a turbulence  trip  wire  called  for 
a wire  diameter  of  about  0.011  inch  (0.28  mm).  This  was  an  odd  size 
and  since  there  was  a large  quantity  of  0.012  inch  (0.305  mm)  piano 
wire  available,  the  latter  was  the  best  choice.  Lengths  of  the  wires 
were  fastened  with  screws  at  the  foil  tips,  and  each  one  was  held  down 
by  counter  sunk  screws  at  two  points  one-fourth  of  the  span  inboard 
from  each  tip.  With  these  few  connection  points,  each  wire  was  free 
to  strum.  This  appeared  to  help  its  effectiveness. 

BLOCK  GAUGE  DYNAMOMETER  AND  INSTRUMENTATION 

It  was  decided  to  measure  the  total  forces  and  moment  directly, 
rather  than  having  to  run  the  test  matrix  twice  in  order  to  determine 
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TABLE  1 - FLAT  PLATE  HYDROFOIL  GEOMETRY 


t - 1.25" 


0.012" 

WIRE  TRIP 


20.125 


.625" 

RADIUS 

312" 


Aspect  ratio 

A - 4 

Thickness  ratio 

t/c  - 0.05208 

Span 

b • 8 ft  * 

2.438  m 

Chord 

c = 2 ft  * 

0.6096  m 

Planform  area 

S - 16  ft^ 

- 1.486  m^ 

Wetted  area  (excluding  portion 
covered  by  bottom  connecting  plate) 

31.68  ft^ 

“ 2.943  m^ 

Ratio  of  foil  wetted  area  to 
planform  area 

1.98 

Cross  section  area  moment  of 
inertia  about  centerline 

1 - 3.468 

inch^  = 1.44  X 10  ^ m^ 

Turbulence  stimulator  wires 

diameter  •» 

0.012  inch  “ 0.305  mm 
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the  strut  tare  forces.  This  was  accomplished  by  positioning  the 
dynamometer  between  the  foil  and  the  base  of  the  supporting  strut. 

A system  of  four,  standard  4-lnch  DTNSRDC  block  gauges  was  arranged 
such  that  two  1000-pound  gauges  and  two  200-pound  gauges  were  used 
for  the  Y-force  and  X-force  measurements,  respectively.  To  help 
reduce  the  influence  of  unmeasured  moments  in  this  system,  pin  joint 
hinges  separated  each  of  the  pairs  of  stacked  X-  and  Y-force  gauges. 

The  hinges  were  located  chordwise  so  that  the  line  of  force  of  each 
of  the  Y-force  gauges  was  equidistant  from  the  quarter-chord  of  the 
foil.  Since  the  expected  center  of  pressure  of  a thin,  uncambered, 
and  unswept  foil  is  at  the  quarter-chord,  each  Y-force  gauge  could  be 
expected  to  share  about  half  the  vertical  force.  The  entire  assemblage 
was  bounded  below  by  a 1.25  inch  (0.03175  m)  thick,  5.25  inch  (0.1334  m) 
wide  bottom  connecting  plate  that  was  bolted  firmly  to  the  hydrofoil  at 
midspan,  and  was  bounded  above  by  a similar  1.25  inch  (0.03175  m)  thick 
upper  connecting  plate  that  in  turn  was  attached  to  the  base  of  the 
standard  9-foot  (2.743  m)  long  strut.  The  robust  nature  of  the 
connecting  plates  gave  the  entire  block  gauge  system  a virtual  stiffness 
that  exceeded  the  stiffness  properties  of  the  individual  block  gauges. 
Figure  2 is  a schematic  of  the  main  parts  of  the  assembled  block  gauge 
dynamometer . 

Instrumentation  for  the  collection  of  the  four  channels  of  force 
data  consisted  of  ENDEVCO  signal  conditioners  for  the  differential 
reluctance  block  gauges  and  DANA  amplifiers.  The  voltage  output  from 
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the  amplifiers  was  digitized  using  ANALOGIC  analog-to-dlgltal  converters. 
All  the  data  were  processed  on-line  using  an  Interdata  Model  4 computer, 
with  automatic  data  storage  on  magnetic  tape  using  a Tri-Data  Cartiflle 
recorder.  Hard  copies  of  the  data  were  available  from  a teletype 
> machine  which  was  used  interactively  to  carry  out  the  data  collection 

procedure. 

SHROUD 

A sheet  metal  shroud  fabricated  from  0.25  inch  (6.35  mm)  thick 
aluminum  enclosed  the  block  gauge  assemblage  and  was  bolted  to  the 
upper  connecting  plate  leaving  a 1/8  inch  (3.175  mm)  gap  between  the 
bottom  edge  of  the  shroud  and  the  upper  surface  of  the  foil.  There 
was  a similar  size  gap  between  the  inner  sides  of  the  shroud  and  the 
side  edges  of  the  bottom  connecting  plate.  The  shroud  had  the  external 
appearance  both  from  the  side  and  from  the  front  of  a rectangular  pod, 
6.125  inches  (0.1556  m)  wide,  10.5  inches  (0.2667  m)  high  including 
the  gap  at  the  bottom  edge,  and  42.125  Inches  (1.07  m)  long  overall. 
Although  this  pod-like  structure  was  not  streamlined  in  the  usual 
sense  of  a nacelle  shape,  it  was  faired  in  the  longitudinal  direction 
with  a 2:1  elliptic  nose  and  a straight  wedge  tail  having  an  included 
angle  of  26  degrees.  Figure  3 is  a sketch  of  the  shroud  geometry. 

Since  the  shroud  was  attached  to  the  upper  connecting  plate  only, 
and  did  not  touch  the  foil  or  bottom  connecting  plate,  no  hydrodynamic 
forces  on  it  were  transferred  to  the  foil  or  to  the  floating  portion  of 


I 
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the  dynamometer.  Thus,  aside  from  hydrodynamic  Interference  effects, 
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the  hydrofoil  was  isolated  from  the  other  parts  of  the  strut-shroud 
system  exposed  to  the  streamflow.  The  volume  inside  the  shroud  was 
flooded  so  the  block  gauges  were  operated  wet.  The  block  gauge  core 
assemblies  and  electrical  cables  were  completely  waterproof,  the  latter 
were  collected  together  in  a bundle  and  led  out  through  the  interior 
of  the  strut  via  access  holes. 

The  shroud  represents  a discontinuity  in  the  strut  geometry. 

For  the  particular  case  of  foil  submergence  f/c  = 0.25  it  pierces 
through  the  free  surface,  and  therefore  plays  the  role  of  the  strut 
in  that  situation. 

STRUT 

The  strut  had  a chord  length  of  30  Inches  (0.762  m) , was 
3 inches  (0.0762  m)  thick,  and  had  a simple  ogive  thickness  variation 
in  the  longitudinal  direction.  The  top  end  of  the  strut  was  attached 
to  a turntable  assembly  that  was  held  in  a support  frame  that  could 
be  raised  and  lowered  on  the  vertical  rails  of  the  towing  carriage. 
Tilting  the  turntable  assembly  provided  the  means  of  changing  angle 
of  attack  of  the  foil.  A portable  inclinometer  was  used  to  measure 
the  angle  of  attack  and  a scale  fixed  to  the  vertical  rails  made  it 
possible  to  fix  the  foil  depth  of  submergence.  Figure  4 is  a sketch 
of  the  entire  foil-strut-support  system  pictured  with  the  hydrofoil 
at  a positive  angle  of  attack.  Note  that  the  Y-force  and  X-force 
values  were  measured  in  a coordinate  system  that  rotated  with  the  foil. 
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CALIBRATION 


The  foil-block  gauge-strut  arrangement  was  hung  on  a test  stand 
for  calibration.  Each  of  the  block  gauges  had  been  calibrated 
individually  before  assembly  so  that  the  sensitivities  of  the  two 
X-gauges  and  the  two  Y-gauges,  respectively,  had  been  preset  as 
closely  as  possible.  The  final  force-to-voltage  slopes  of  the 
complete  dynamometer  were  determined,  however,  in  the  assembled 
condition. 

ERRORS  AND  INTERACTIONS 

Inaccuracies  in  the  force  and  moment  data  are  inherent  to  the 
mechanical/electrical  system  employed  and  are  also  caused  by  load 
distortions  of  both  the  model  and/or  force  gauge  system.  Quantitative 

i 

estimates  of  these  errors  are  outlined  here. 

Measurement  errors  in  the  values  of  the  forces  manifest  themselves 
in  small  scatter  in  the  instrument  readings  of  the  block  gauge  outputs 
during  calibration.  For  both  the  X-  and  Y-forces,  these  errors 
appeared  to  be  always  less  than  + 0.01  of  the  applied  load  over  the 

entire  load  range.  Angle  of  attack  values  a were  determined  to  | 

j 

within  + 0.02  degrees  so  that  the  relative  error  at  a>2  degrees  was  I 

less  than  or  equal  to  + 0.01.  Depth  of  submergence  was  accurate  to 

within  about  +0.01  foot  (+  0.003048  m)  so  that  for  the  submergence  | 

range  tested,  the  relative  angle  error  was  less  than  or  equal  to  + 0.04. 

Elastic  deflections  of  the  model  hydrofoil  under  steady  hydrodynamic 
loads  undoubtedly  altered  the  geometry  of  the  lifting  surface.  Two 
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categories  of  simple  distortions  have  been  estimated:  vertical 
deflection  associated  with  the  spanwlse  bending  caused  by  the 


Y-force  loading,  and  twisting  deflection  due  to  the  hydrodynamic 
pitching  moment.  The  case  of  bending  deflection  has  been  computed 
using  the  assumption  that  an  elliptic  distribution  of  vertical 
loading  acted  over  the  entire  hydrofoil  span  and  that  the  hydrofoil 
half-span  was  a straight  beam  of  constant  section  properties, 
cantilevered  from  the  edge  of  the  bottom  connecting  plate.  The 
maximum  deflection  that  occurred  at  each  wing  tip  under  the  most 
extreme  loading  encountered  was  predicted  using  elementary  strength 
of  materials  to  be  about  0.424  inch  (0.01077  m)  which  is  approximately 
1/3  the  plate  thickness. 

The  torsional  deflection  was  estimated  using  an  approximate 

method  for  computing  the  angle  of  twist  for  rectangular  cross  section 

9 

bars,  as  outlined  by  Seely  and  Smith.  The  pitching  moment  taken 
about  the  quarter-chord  axis  was  transferred  to  the  mid-chord  axis, 

; and  was  assumed  to  be  distributed  elliptically  across  the  span.  For 

A 

I a rectangular  beam  of  constant  section  properties  having  the  same 

I 

i thickness-to-chord  ratio  as  the  hydrofoil  model,  the  most  extreme 

' test  condition  was  predicted  to  produce  a maximum  wing  tip  twist  of 

I 

about  0.24  degrees. 

Elastic  distortions  within  the  force  gauge  arrangement  itself 
were  also  responsible  for  errors  in  the  final  measured  results,  since 
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each  block  gauge  is  in  fact  a deflection  measuring  device.  Unlike 
the  spurious  character  of  the  calibration  scatter  mentioned  previously, 


the  interaction  errors  of  the  block  gauge  were  uni-directional.  That  I 

is,  the  relative  errors  were  always  plus  or  always  minus,  rather  than  | 

the  + values  typical  of  calibraticn  error  bounds.  Interactions  of 

Y-force  into  the  X-readings,  and  X-forces  into  the  Y-readlngs  were  | 

I 

carefully  determined  during  the  calibration  phase  of  the  test.  The  | 

Influence  of  Y-force  on  the  X-readlngs  was  found  to  be  the  larger, 

so  the  strongest  interaction  errors  occurred  for  the  foil  at  the 

largest  angles  of  attack.  Typical  maximum  values  of  the  net  relative 

corrections  predicted  for  the  most  extreme  test  conditions  were  about 

-0.006  for  the  lift  force  and  +0.015  for  the  drag  force.  Since  the  | 

Interaction  slopes  were  themselves  subject  to  calibration  errors  | 

roughly  of  the  order  of  + 0.01  of  the  applied  load  (same  as  the  pure  | 

force  calibration  errors),  it  was  felt  that  such  small  adjustments 

to  the  measured  values  were  not  justified  at  this  time.  ' 

None  of  the  data  presented  in  this  report  have  been  corrected 
or  altered  for  any  of  the  errors  described  above. 

EXPERIMENTAL  PROCEDURES 

Experiments  were  carried  out  on  Carriage  II  in  the  David  Taylor 
Model  Basin  at  the  David  W.  Taylor  Naval  Ship  R&D  Center,  Carderock 
facility.  This  deep  water  towing  tank  has  a water  depth  of  22  feet 
(6.706  m)  with  a width  of  51  feet  (15.54  m) . The  measured  water 
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temperature  was  a constant  68“F  (20'’C),  so  the  density  and  kinematic 

3 

viscosity  values  used  in  data  reduction  were  p = 1.9367  slug/foot 
(998.13  kg/m^)  and  v = 1.0836  x lO"^  foot^/second  (1.0067  x 10~^  m^/s), 


respectively.  The  nominal  test  matrix  for  the  data  presented  in  this 
report  is  given  in  Table  2.  At  the  low  speeds  of  these  experiments, 
no  cavitation  was  ever  observed  to  occur,  even  at  the  highest  angles 
of  attack. 

A typical  test  run  was  made  with  fixed  submergence  depth  f and 
angle  of  attack  a.  Several  passes  down  the  basin  were  necessary  to 
complete  all  the  speeds  desired.  After  zeros  were  collected  at  the 
beginning  of  each  pass,  the  carriage  was  brought  up  to  speed  and 
five  channels  of  data  were  taken;  the  velocity,  and  the  output 
voltages  of  Y^,  and  Y^.  These  were  collected,  averaged  over 

a 10  second  continuous  record  for  each  measured  point,  and  processed  | 

onboard  the  carriage  for  both  teletype  output  and  tape  storage  | 

as  described  earlier.  i 

) 

When  changing  angle  of  attack,  appropriate  changes  had  to  be  ^ 

made  for  the  vertical  setting  of  the  support  bracket  along  the  vertical 
rails  so  as  to  maintain  a constant  depth  of  submergence  from  the 
undisturbed  free  surface  to  the  foil  quarter-chord. 

Since  the  X-  and  Y-forces  were  measured  in  the  foil  coordinate 
system,  they  were  resolved  into  the  desired  drag  and  lift  components 
(forces  parallel  and  perpendicular  to  the  free  stream  direction) 
using  the  formulas 
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TABLE  2 

NOMINAL  TEST  MATRIX  FOR  LOW  FROUDE  NUMBER 
HYDROFOIL  EXPERIMENTS 


I Parameter  Range 


1.  Speed 

U » 4,  6,  8,  10,  12,  16, 

20,  24,  28,  ft/sec 

U = 1.22,  1,83,  2.44,  3.05, 

3.66,  6.1,  7.32,  8.53,  m/s 

, Chord  Froude  Number 

F = 0.5,  0.75,  1.0,  1.5, 

^ 2.0,  2.5,  3.0,  3.5 

Chord  Reynolds  Number 

R /lO^  = 0.738,  1.11,  1.48, 
1.85,  2.21,  2.95, 
4.43,  5,17 

2.  Angle  of  Attack 

a = 0,  2,  4,  6,  8 degrees 

3.  Depth  of  Submergence 

(to  foil  quarter-chord) 

f = 0.5,  1.0,  1,5,  2.0,  3.0, 
4.0,  7.0  ft 

f = 0,152,  0.305,  0.457,  0.61, 
0.91,  1.22,  2.13  m 

Snbmereence  Ratio 

f/c  = 0.25,  0.5,  0.75,  1.0, 
1.5,  2.0,  3.5 
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D - (Xj^  + X2)cos  a + (Y^  + Y2)sln  a 

L = (Yi  + Y2)cos  a - (X^  + X2)sln  a 


(1) 


where  here  (X^  + X^)  and  (Y^  + Y^)  denote  the  sums  of  the  block 
gauge  outputs  converted  to  forces  using  the  calibration  slopes. 

The  hydrodynamic  pitching  moment.  In  foot-pounds,  measured  about 
the  foil  quarter-chord  (positive  for  nose  up)  was  determined 
from  the  expression 

M(c/4)  " 0.708333  (Y^  - Y^)  (2) 


where  Y^  and  Y^  represent  the  weighted  force  values  of  the 
forward  and  after  Y-force  gauges,  respectively.  The  two  force- 
to-voltage  slopes  for  Y^  and  Y^  were  determined  during  calibration 
to  satisfy  simultaneously  both  the  applied  moment  as  well  as  the 
net  applied  vertical  force,  and  therefore  properly  account  for  slight 
differences  in  the  actual  sensitivities  of  the  two  Y-gauges  in  the 
assembled  system. 

The  distance  to  the  center  pressure  x , measured  from  the 

cp 

leading  edge  of  the  foil  was  determined  from 


X 

cp 


(3) 


Nondlmensional  coefficient  forms  of  the  drag,  lift,  moment 
and  center  of  pressure  location  are  defined  in  the  usual  fashion  as 
noted  in  the  list  of  notation. 
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RESULTS  AND  DISCUSSION 


Tables  of  the  measured  force  and  moment  data  reduced  to 
coefficient  form  are  presented  in  Appendix  A.  For  purposes  of 
this  discussion,  a considerable  amount  of  cross  plotting  and 
further  analysis  has  been  carried  out. 

EFFECT  OF  FROUDE  NUMBER  ON  DRAG  AND  LIFT 

Curves  of  total  measured  drag  coefficient  and  lift  coefficient 
are  plotted  versus  chord  Froude  number  in  Figures  5 through  11  for 
constant  submergence  ratios  of  f/c  = 0.25  through  3.5,  respectively. 

The  contours  are  for  constant  geometric  angle  of  attack  a.  It  is 
evident  that  as  the  foil  operates  closer  to  the  surface,  the  effects 
of  speed  (Froude  number)  become  more  and  more  pronounced.  The  curves 
of  C versus  Froude  number  are  nearly  flat  at  the  deeper  submergences 
(as  they  should  be).  For  the  foil  operated  nearer  the  surface,  the 
dip  in  lift  coefficient  in  the  vicinity  of  F^  ~ 1 becomes  increasingly 
exaggerated . 

At  the  shallowest  submergence  of  f/c  = 0.25,  the  hydrofoil 
experiences  negative  lift  throughout  the  tested  speed  range  at  a 
geometric  angle  of  attack,  a = 0,  and  displays  intervals  of  negative 
lift  for  certain  Froude  numbers  even  up  to  an  angle  of  attack  of 
a = A°.  For  the  a = 0 case,  the  miminum  point  of  the  C curve  has  a 

Li 

value  of  C,  = -0.21  at  a Froude  number  F =0.75  (F,  = 1.5).  With 
L c t 

the  f/c  = 3.5  curves  of  Figure  11  as  a measure  of  wave-free  performance. 
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It  can  be  seen  that  this  foil  requires  nearly  a 4 degree  angle  of 


attack  In  order  to  reach  a positive  = 0.21.  Therefore,  at  f/c  = 

0.25  for  a speed  near  F * 0.75,  the  foil  apparently  sees  a net  or 

c 

averaged  flow  angle  of  almost  minus  4 degrees.  It  is  not  known 
definitely  by  separate  experiment  or  by  calculation  whether  such 
a downwash  flow  could  be  generated  by  a 5.2  percent  thick  plate  alone. 

At  the  shallowest  submergence  of  f/'c  = 0.25,  the  plate  thickness-to- 
depth  ratio  is  only  t/f  = 0.2083,  and  the  foil  probably  experiences 
some  wavemaking  drag  due  to  thickness  as  well  as  wave-induced 
negative  velocities  produced  by  thickness  across  the  entire  span. 

Of  course,  strut-  and  pod-induced  downwash  must  alter  the  flow 
angle  of  attack  locally  near  the  center  span  region.  However,  it  seems 
unlikely  that  a net  minus  4 degrees  could  be  induced  across  the  8 foot 
(2.44  m)  span.  Unfortunately,  there  are  no  computation  procedures 
currently  available  for  estimating  the  strut-  and  pod-induced  downwash 
velocities  at  the  intermediate  values  of  Froude  number  of  interest  here. 
Such  calculations  could  allow  the  best  possibility  of  isolating  the 
foil-alone  performance  from  the  present  data. 

The  stong  dip  in  the  variation  with  speed  for  the  smallest 
submergence  cases  is  present  for  all  values  of  a tested,  with  the 
curves  moving  up  in  parallel  contours  for  various  angles. 

Regarding  the  drag  coefficient  variation  versus  speed  in 
Figures  5 through  11,  it  can  be  seen  that  for  operation  nearer  and 
nearer  the  surface  the  drag  coefficient  curves  steepen  considerably. 
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This  reflects  the  increasing  amount  of  hydrofoil  wavemaking  drag. 

It  is  reasonable  to  suspect  that  a portion  of  this  is  caused  by 
wavemaking  due  to  thickness  effects  and  some  interference  effects 
due  to  induced  downwash  flow  from  the  strut  and  pod. 

For  reference,  in  each  of  Figures  5 through  11,  the  1957 
International  Towing  Tank  Conference^^  (ITTC)  turbulent  friction 
correlation  line  is  plotted,  converted  to  a friction  drag  coefficient 
based  on  the  planform  area. 

EFFECT  OF  FOIL  SUBMERGENCE  ON  DRAG  AND  LIFT 

Curves  of  total  measured  drag  and  lift  coefficients  are  plotted 

versus  submergence  ratio  in  Figures  12  through  15  for  the  constant 

chord  Froude  numbers  of  F = 0.5,  0.75,  1.0,  and  2.0,  respectively. 

c 

The  contours  are  for  constant  angle  of  attack.  Interesting  changes 

In  the  characteristics  occur  at  the  low  Froude  numbers.  These  are 

especially  noticeable  in  comparing  the  F^  = 0.5  case  with  the  other 

three  cases.  The  F =2  curves  in  Figure  15  are  typical  of  the 
c 

expected  behavior  of  drag  and  lift  coefficients  with  depth  of 

submergence  at  large  Froude  numbers:  both  force  coefficients  fall 

off  with  decreasing  submergence  depth,  with  the  lift  coefficient 

suffering  the  most  rapid  decrease.  As  the  Froude  number  is  decreased 

down  to  F *0.75,  the  characteristic  curves  are  similar  to  the  F =2 
c c 

case,  but  display  tendencies  to  remain  level  or  have  slight  humps  in 
the  drag  and  lift  coefficients  before  the  final  rapid  descent  at 
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0.5  there  is  a definite 


submergence  ratios  below  0.75.  At  F ■ 

c 

difference.  The  drag  coefficients  actually  increase  with  small 

f/c  ratios,  and  the  lift  coefficients  rise  to  rather  sharp  peaks 

before  falling  off  below  f/c  < 0.5.  Evidently  an  important  change 

in  the  flow  occurs  for  shallow  submergences  in  the  Froude  number 

range  between  F = 0.5  aiid  0.75.  From  visual  observations  of  the  free 
c 

surface,  this  change  is  manifest  in  a spanwlse  breaking  wave  or 
hydraulic  jump  that  appears  over  the  foil  in  this  speed  range.  The 
photograph  Included  in  the  frontispiece  of  this  report  shows  the 
spanwlse  wave  disturbance  observed  at  f/c  = 0.25  and  F^  = 0.5.  The 
dip  the  C versus  F curves  mentioned  earlier  also  occurs  in  the 

L C 

same  Froude  number  range  as  the  appearance  of  the  breaking  waves. 

It  is  Interesting  to  note  at  this  point  that  in  two-dimensional 

hydrofoil  experiments  described  by  Parkin,  et  al,^^  similar  waves 

were  encountered  and  were  shown  to  have  a drastic  effect  on  the 

pressure  distribution  over  the  foil  surface.  For  the  particular  case 

of  a 12  percent  thick  Joukowskl  profile  hydrofoil  at  a = 5 degrees 

and  with  the  trailing  edge  submerged  a distance  of  0.25c,  the 

character  of  the  measured  pressure  distributions  In  Reference  11  showed 

a remarkable  transformation  at  about  F =0.61.  For  Froude  numbers 

c 

higher  than  0.61  the  pressure  distributions  resembled  those  of  a 
hydrofoil  at  deep  submergence.  At  Froude  numbers  below  F^  = 0.61,  the 
upper  surface  negative  pressure  peak  was  suddenly  shifted  to  near  the 
rear  of  the  foil  with  a strong  positive  pressure  occurring  at  the  nose. 
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This  dramatic  change  was  definitely  associated  with  a hydraulic  jump 
that  appeared  over  the  hydrofoil. 


L 


1 

I 


In  the  present  finite  aspect  ratio  experiments,  the  apparent 

shift  in  flow  character  occurs  between  F “0.5  and  0.75  and  the 

c 

effects  are  noticeable  In  the  character  of  the  net  hydrodynamic 
forces.  While  Interference  effects  from  the  strut  and  pod  could 
account  for  some  of  the  changes.  It  appears  that  the  primary 
Influence  Is  the  breaking  wave  surface  disturbance  which  was  observed 
to  occur  across  the  entire  span  of  the  foil. 

EFFECT  OF  ANGLE  OF  ATTACK  ON  DRAG  AND  LIFT 

Graphs  of  drag  and  lift  coefficients  plotted  versus  angle  of 

attack  are  given  in  Figures  16  through  21  for  constant  submergence 

ratios  of  f/c  = 0.25,  0.5,  0.75,  1.0,  1.5,  and  2.0,  respectively.  The 

contours  are  for  selected  values  of  chord  Froude  number  F = 0.5,  0.75, 

c 

1.0,  and  3.0.  Froude  number  effects  are  evident  at  the  low  speeds  in 
the  separation  of  curves  and  variation  of  lift-curve  slope  for  various 
Froude  numbers. 


MOMENT  AND  CENTER  OF  PRESSURE 


Curves  of  the  pitching  moment  coefficient  about  the  quarter-chord 
Cj^  and  the  center  of  pressure  ratio  are  plotted  versus  chord 

Froude  number  in  Figures  22  through  28  for  constant  submergence  ratios 
of  f/c  = 0.25  through  3.5,  respectively.  The  contours  are  for  constant 
geometric  angle  of  attack. 
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r- 1 


As  expected,  the  Influence  of  the  Froude  number  is  more  dramatic 

for  the  shallow  submergence  cases,  with  the  moment  coefficient  having  a 

strong  dip  centered  near  F •=  0.5  - 0.75.  The  minimum  point  moves  to 

c 

higher  Froude  numbers  with  increasing  f/c.  The  noteworthy  changes 

in  the  sign  and  magnitude  of  the  moment  in  the  Froude  number  range 

of  F = 0.5  - 1.0  are  associated  with  the  observed  spanwide  breaking 
c 

wave,  and  can  be  explained  qualitatively  in  terms  of  the  pressure 
distribution  results  mentioned  earlier. 

The  center  of  pressure  location  also  reflects  the  exaggerated 
Froude  number  effects  in  the  vicinity  of  F^  = 0.5  to  1.5.  Movement 
of  the  center  of  pressure  off  the  foil  is  a common  feature  of  the 

12 

effect  of  camber  at  low  angles  of  attack  (see,  for  example  Von  Mlses  ). 

Although  the  foil  itself  is  uncambered,  there  are  clearly  vertical 

velocities  induced  by  the  free  surface  and  by  the  presence  of  the 

strut  and  pod.  These  create  a Froude-dependent  camber-like  effect. 

For  the  deepest  submergence  case,  in  Figure  28,  the  center  of 

pressure  ratio  appears  to  settle  down  roughly  in  the  neighborhood  of 

the  anticipated  value  of  x tc  - 0.25,  although  there  is  still  an 

cp 

effect  of  angle  of  attack  that  moves  the  center  of  pressure  aft  for 

increasing  a.  It  is  Interesting  that  for  the  larger  values  of  a,  the 

curves  of  center  of  pressure  ratio  versus  Froude  number  seem  to  level 

off  at  roughly  the  same  value  of  x /c  = 0.22  even  at  the  shallower 

cp 

submergences,  and  despite  the  dip  occurring  near  F = 0.75. 

c 
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LIFT-CURVE  SLOPE 


Figures  29  and  30  show,  respectively,  the  Influence  of  Froude 

number  and  depth  of  submergence  on  the  lift-curve  slope.  The  curves 

of  C,  versus  F with  contours  of  submergence  ratio  in  Figure  29 
L c 

a 

indicate  that  in  the  low  Froude  number  range  near  F ■ 1,  the 

c 

corresponding  dips  in  the  plots  of  C versus  F are  brought  about 

L C 

not  only  by  an  overall  shift  in  the  net  angle  of  attack  seen  by  the 
foil,  but  also  by  a reduction  in  the  lift-curve  slope.  This  is 
probably  caused  by  a camber  or  flow  curvature  effect  due  to  the 
observed  spanwise  broken  waves  occurring  in  this  speed  range  at  the 
shallowest  submergences,  and  due  to  Interference  velocities  from  the 
strut  and  pod.  In  any  case,  this  graph  shows  that  in  addition  to 
submergence  ratio,  Froude  number  must  also  be  considered  in  determining 
lift-curve  slopes. 

In  Figure  30,  the  reference  value  of  lift-curve  slope  for  the 
present  data  is  that  for  the  deepest  submergence  and  highest  speed 

C = C (f/c  = 3.5,  F = 3.5)  *=  0.0617  deg  (4) 

n(REF)  a 

The  solid  line  curve  in  Figure  30  is  taken  from  the  experimental 

results  of  Wadlin,  et  al^  for  a rectangular  planform,  aspect  ratio  4 

hydrofoil  at  chord  Froude  numbers  F^  > 3.24.  The  distorting  effect 

of  small  Froude  number  is  evident  here  in  the  curves  for  F * 1 and  2. 

c 

There  Is  good  correspondence  between  the  present  lift-slope  ratio 

results  at  F » 3 and  the  results  of  the  previous  experiments, 
c 


COMPARISON  OF  LIFT  RESULTS 


The  low  speed  end  of  the  NACA  experiments  of  Reference  5 can 

also  be  used  in  a comparison  of  lift  at  similar  Froude  numbers  and 

submergences.  Comparison  plots  of  lift  coefficient  versus  angle  of 

attack  measured  from  zero  lift  are  given  in  Figures  31  through  34 

with  the  present  data  at  F • 3.0  and  f/c  * 0.5,  1,  2,  and  3.5, 

c 

respectively.  The  NACA  data  are  for  a constant  F^  = 3.24  and 
f/c  = 0.59,  1.09,  2.09,  and  3.09/4.09,  respectively.  It  can  be  seen 
that  absolute  values  of  lift-curve  slopes  for  the  flat  plate  foil 
appear  to  be  somewhat  less  than  those  for  the  NACA  64A412  section 
used  in  Reference  5.  The  discrepancies  may  be  due  to  the  poor  section 
properties  of  the  faired  flat  plate,  but  the  differences  in  the  C^ 

a 

values  decrease  with  increasing  submergence  to  about  6 percent  at 
the  deepest  cases  shown. 

RESIDUAL  DRAG  COEFFICIENTS 

An  attempt  has  been  made  to  extract  the  induced-plus-wave  drag 
from  the  measured  values  of  total  drag.  In  the  present  work,  the 
residual  drag  coefficient  is  defined  as 


(R  ) 

R D D . c 
vise 


(5) 


where  C is  the  estimated  total  viscous  drag  given  by 

vise 


=1)  , "c*  - \ «c>  + S 

vise  f vp 


(6) 


i 
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Here  C denotes  the  flat  plate  friction  drag  coefficient  based  on 
°f 

planfonn  area,  R Is  the  foil  chord  Reynolds  number,  and  C_  denotes 
c IJ 

vp 

the  viscous  pressure  or  "fonn"  drag.  It  was  decided  to  use  the  ITTC 

1957  friction  correlation  curve  for  the  flat  plate  friction  drag 

Instead  of  the  usual  practice  of  using  the  Schoenherr  line.  Recent 
13 

work  by  Gpanvllle  has  shown  that  at  low  turbulent  flow  Reynolds 

numbers  (5  x 10^  < R < 10^),  the  ITTC  line  turns  out  to  be  very  close 

c 

to  the  best  semi- empirical  flat  plate  friction  line.  It  is  therefore 
preferable  to  the  Schoenherr  curve  in  the  range  of  the  present 
experiments  and  becomes  indistinguishable  from  the  Schoenherr  values 
at  Reynolds  numbers  greater  than  10^. 

For  the  deepest  submergences  tested  (f/c  = 2 and  3.5),  at  the 
highest  speeds,  and  at  zero  angle  of  attack  the  drag  coefficients  are 
observed  to  level  out  on  a curve  parallel  to  the  ITTC  1957  line.  Since 
at  these  depths  the  C curve  for  a •»  0 is  flat  and  very  near  zero,  the 

JL 

difference  between  the  measured  total  drag  coefficient  and  the  flat 
plate  friction  value  is  assumed  to  be  the  viscous  pressure  drag, 
entirely  free  of  lift-dependent  wavemaking  effects.  The  result  is 
an  estimate  for  the  form  drag  coefficient  of 


Cp  = 0.005 
vp 

This  is  assumed  to  be  a function  of  shape  only,  independent  of 

Reynolds  number  so  that  all  the  Reynolds  number  variation  in 

is  contained  in  C of  Equation  (6). 

^f 


(7) 


vise 
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Inclusion  of  the  viscous  pressure  drag  in  the  definition  of 

residual  drag  is  somewhat  unusual,  since  the  typical  practice  is 

to  absorb  any  small  form  drag  into  the  total  pressure  drag  category. 

The  definition  in  Equation  (6)  is  desirable  here  because  the  idea  of 

the  present  experiments  was  to  remove  the  section-shape-dependent 

drag  as  completely  as  possible,  leaving  only  the  Froude-dependent  drag 

variation.  The  form  drag  of  this  section  was  determined  to  be  rather 

large,  and  therefore  an  important  contributor  to  total  section  drag. 

This  is  not  an  unexpected  result  for  the  drag  on  a faired  flat  plate 

of  the  type  used  here.  One  set  of  examples  can  be  found  in  a NACA 

14 

report  by  Wadlln,  et  al,  with  results  of  experiments  on  three  flat 
plate  hydrofoils  of  small  aspect  ratio  having  2:1  elliptical  noses 
and  straight  wedge  tails,  with  thickness  ratios  of  1.0,  1.3,  and 
2.7  percent.  The  measured  drag  coefficients  at  the  deepest  submergence 
cases,  taken  at  zero  lift,  were  typically  0.004  to  0.0025  above  the 


ITTC  1957  friction  line  in  the  same  Reynolds  number  range  of  the  present 
experiments.  Such  values  of  form  drag  for  very  thin  plates  lend 


confidence  in  the  value  determined  in  the  present  experiments. 

vp  2 

Curves  of  measured  C /C  versus  chord  Froude  number  are  plotted 

R L 

in  Figures  35  through  41  for  submergence  ratios  of  f/c  = 0.25  through 


3.5,  respectively.  The  contours  are  for  constant  angle  of  attack.  As 


noted  earlier,  the  residual  drag  coefficients  C should  be  a close 

R 

approximation  to  the  Froude-dependent  variation  of  (C  + C,,)  which 

^1  ” 

is  the  sum  of  total  induced  drag  plus  wavemaking  drag.  For  a hydrofoil 
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1 


near  the  free  surface,  the  total  Induced  drag  consists  of  the 
unbounded  flow  plus  the  biplane  induced  drag. 

2 

For  reference,  the  proper  limits  of  + C^)/C^  at  Infinitely 
large  and  zero  Froude  numbers  are  indicated  in  the  Figures  35  through 
41  as  horizontal  dash-dot  lines  at  the  extremes  of  the  speed  range. 

In  the  limit  of  infinitely  large  Froude  number,  the  wavemaking  drag 
disappears  and  all  that  remains  is  the  total  Induced  drag  including 
the  biplane  effect  of  the  free  surface.  Thus 


2 


1 + 5 + o(X) 
ttA 


(8) 


where  6 = Glauert  planform  factor  for  non-elliptic  planforms^^ 

a(X)  = Prandtl  biplane  factor,  which  is  strictly  a function 
of  the  submergence-to-semi  span  ratio  X 

A discussion  of  the  biplane  function  a(X)  can  be  found  in  von  Karmen 

and  Burgers. 

In  the  limit  of  zero  Froude  ntimber,  the  free  surface  acts  as  a 
rigid  wall.  The  sign  of  the  biplane  Induced  drag  must  be  reversed 
from  the  high  Froude  number  case,  since  the  flow  is  now  s)mimetrlcal 
about  the  free  surface  plane.  For  vanishing  Froude  number,  the 
wavemaking  drag  goes  to  zero,  leaving  the  total  induced  drag  at  the 
low  speed  extreme  expressed  as 


2 


1 + 5 - g(X) 
ffA 


(9) 
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The  factor  6 Is  equal  to  0.03  for  a rectangular  planform,  aspect  ratio  4 


lifting  foil.  Numerical  values  of  o(X)  valid  for  an  elliptic  load 

2 

distribution,  and  values  of  the  two  limits  of  (C„  + C,,)/C,  that 

w L 

appear  in  Figures  35  through  41  are  given  in  the  table  of  Appendix  B. 

A rough  comparison  between  a theoretical  prediction  of 
2 2 

(C  + C )/C  and  measured  C /C  versus  F is  shown  in  Figure  42. 

U , W l_i  tv  i_i  C 

^ 4 

The  theory  curve  is  from  Nishlyama  for  A = 4,  and  applies  to  a sub- 
mergence ratio  of  f/c  = 0.59.  The  two  experimental  curves  are  for 
f/c  = 0.5  and  0.75,  at  angle  of  attack  a = 8°  only.  This  indicates 
that  at  least  for  large  lift,  there  is  a measure  of  correspondence 
between  the  experiments  and  the  analytical  results  of  linear  potential 
theory.  The  same  cannot  be  said  for  smaller  angles  of  attack. 

It  can  be  seen  from  Figures  35  through  41  that  the  curves  of 
2 

C„/C,  versus  F for  various  angles  of  attack  do  not  fall  on  the 
R L c 

same  curve,  although  the  spread  becomes  considerably  smaller  as  the 
depth  of  submergence  is  increased,  and  at  the  larger  angle  of  attack 

for  any  given  submergence.  The  prediction  of  linear  potential  theory 

3 4 2 

(Breslin  , Nishlyama  ) is  that  there  is  one  curve  of  (C„  + C,,)/C, 

w L 

for  an  Isolated  hydrofoil  of  a given  planform  shape  and  depth  of 
submergence,  regardless  of  angle  of  attack.  One  explanation  of  the 
divergence  between  the  experimental  results  and  the  theory  is  that 
the  linear  theory  may  be  Inadequate  to  describe  the  flow  phenomena 
at  low  chord  Froude  numbers,  especially  at  the  shallowest  submergences 
where  spanwide  breaking  waves  are  observed  to  occur  in  just  the  Froude 
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number  range  of  the  most  exaggerated  differences.  Another  possibility 

is  that  the  present  measurements  Include  Froude-  and  depth-dependent 

interference  drag  ''omponents  caused  by  the  strut  and  pod. 

A complementary  viewpoint  of  the  variation  of  residual  drag 

coefficient  is  afforded  in  Figures  43  through  46  for  submergence 

ratios  of  f/c  = 0.25,  0.5,  0.75,  and  1.0,  respectively.  In  these 

graphs,  C is  plotted  versus  positive  lift  coefficient  squared 

to  the  right  of  zero,  and  versus  negative  lift  coefficient  squared 

to  the  left.  The  contours  are  for  chord  Froude  number  equal  to  0.5, 

0.75,  1.0,  and  3.0.  From  the  linear  theory,  the  expected  variation 

2 

of  (C  + C,,)  versus  (+  C, ) for  a given  Froude  number  consists  of 
w — L 

a straight  line  in  each  quadrant,  extending  from  zero  in  the  positive 

C direction.  For  large  C values,  the  curves  do  appear  to  settle 
R L 

onto  straight  lines.  Near  zero  lift,  however,  there  is  a mixture  of 

nonlinear  lift-,  Froude-,  and  depth-dependent  drag  revealed  by  the 

rapidly  changing  slopes  and  the  non-zero  values  of  C at  C =0. 

R L 

These  drag  effects  are  accentuated  at  shallow  submergences,  and  seem 
to  diminish  in  magnitude  as  the  submergence  depth  is  increased. 

Although  it  is  not  possible  to  separate  out  all  the  drag 
components  that  are  present  in  the  C values  of  these  experiments, 
it  is  fair  to  say  that  the  drag  due-to-llft  of  a hydrofoil  at  low 
Froude  numbers  is  distinctly  influenced  by  the  presence  of  a strut 
and  pod.  These  factors  should  be  better  understood  both  from  a point 
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of  view  of  designing  future  experiments  and  for  estimating  drag 
on  large  prototype  hydrofoil  systems. 


APPLICATION  OF  THE  MEASURED  DATA 

An  example  application  of  the  model  drag  and  lift  data  is 

presented  here  for  the  performance  of  a hypothetical  large  hydrofoil 

system  consisting  of  three,  aspect  ratio  4,  rectangular  planforms 

operated  at  constant  lift.  If  a length  scaling  ratio  of  9 Is  used, 

then  the  prototype  chord  length  is  18  feet  (5.486  m) , and  the  prototype 

speed  interval  corresponding  to  a chord  Froude  number  range  of  F = 

c 

1.75  to  3.5  is  Vj,  = 25  t^^O  knots.  This  is  a reasonable  speed  range 
for  operation  between  takeoff  and  cruise. 

Several  questions  can  be  investigated: 

1.  What  is  the  predicted  breakdown  of  drag  for  a large-chord 
hydrofoil  operated  at  constant  lift  in  a practical  range 
of  subcavitat ing  speeds?  Is  the  wavemaking  portion  of 
the  drag  significant? 

2.  How  does  the  trend  of  measured  residual  drag  versus  Froude 
number  compare  with  the  trend  of  induced  drag  alone? 

3.  What  is  the  effect  of  lift  loading  on  the  residual  drag?  ) 

iWlth  the  available  lift  coefficients  of  these  experiments,  the  | 

2 ^ 

attainable  lift  loadings  are  limited  to  about  L/S  = 800  pound/foot  j 

r 2 

(38304  N/m  ) or  smaller,  in  the  speed  range  V < 50  knots.  Unfortunately, 

K 

this  does  not  match  with  the  usual  requirements  of  L/S  = 1200  to  1400  I 

2 2 ' 
pounds/foot  (57456  to  67032  N/m  ),  but  such  lift  loadings  can  only 

be  achieved  with  cambered  foils  and/or  with  flaps.  This  illustrative 

example  is  worked  using  only  the  measured  data. 
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1 


The  case  of  f/c  “ 1 Is  considered  as  representative.  With  a 

2 

lift  loading  of  L/S  « 800  pounds/foot  , the  total  weight  supported  by 

2 

all  three  hydrofoils,  each  one  having  a planform  area  of  1296  foot 
2 

(120.4  m ),  Is  1389  long  tons  (1411  m. tonne).  Measured  model  lift 

coefficients  and  deduced  residual  drag  coefficients  are  assumed  to 

apply  at  the  same  Froude  numbers  for  the  prototype.  Lift  coefficients 

corresponding  to  constant  lift  loading  on  each  foil  are  calculated, 

and  cross  plots  of  C and  C versus  a are  used  to  find  the  required 

L R 

C values  at  constant  lift  by  interpolation.  Friction  drag  coefficients 
at  the  prototype  Reynolds  numbers  are  determined  from  the  ITTC  1957 
correlation  line.  Wavemaking  drag  is  found  by  deducting  the  total 
induced  drag  from  the  residual  drag  as  follows 


where 

r 2 2 

C = L [1  + 5 + o(X)]  = (0.1004)c/  for  A = 4,  f/c  = 1 
°i  irA  ^ 

Then  the  total  projected  drag  coefficient  for  each  hydrofoil  is 

estimated  by 

Table  3 contains  the  numerical  results  of  this  procedure,  with  the 
drag  coefficient  breakdown  versus  Froude  number  plotted  in  Figure  47. 
It  can  be  seen  from  this  extrapolation  of  measured  model  drag  data 
that,  the  wavemaking  drag  is  a significant  portion  of  the  residual 
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drag,  and  of  the  total  hydrofoil  drag  as  well,  especially  at  low 
Froude  numbers. 


r 


The  drag  of  three  foils  Is  taken  as  simply  three  times  the  drag 
on  one.  The  resulting  llft-to-drag  ratio  for  the  foils  alone  Is 
plotted  versus  Froude  number  In  Figure  47.  There  Is  a noticeable 
reduction  of  llft-to-drag  ratio  at  the  lower  speeds  for  operation 
at  constant  lift  at  f/c  = 1. 

It  Is  Interesting  to  compare  the  relative  trends  of  the  residual 
drag  and  the  Induced  drag  as  functions  of  speed.  Figure  48  Is  a 
plot  of  C /C  (F  = 3.5)  and  C /C  (F  = 3.5)  versus  chord  Froude 

i\  O ^ 13  ^ O 

number  for  constant  f/c  = 1.  The  reference  values  for  C and  C are 

R 

taken  at  V = 49.9  knots  (F  = 3.5)  for  a c = 18  ft  hydrofoil.  With  the 

K.  C 

simple  expressions  for  C given  previously,  the  Induced  drag  ratio  at 

°1 


constant  lift  loading. 

for  any  speed  V , 
K. 

is 

(49.9  ^ 

4 

/ 

I 

Cn  (F  = 

3.5) 

V 

F 

°1  " 

1 1 

\ 

(12) 


The  C values  were  obtained  at  various  lift  loadings  by  simple 
R 

Interpolations  as  described  earlier.  The  contours  are  for  prototype 

2 

values  of  L/S  = 600,  800,  and  900  pounds/foot  . Figure  48  Indicates 
that  the  relative  rise  of  measured  residual  drag  is  slightly  milder 
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than  the  relative  rise  of  induced  drag,  at  least  for  the  prototype 

loadings  covered  by  the  present  data.  It  is  also  noted  that  the 

relative  rise  of  the  measured  C steepens  for  larger  lift  loadings. 

R 

CONCLUSIONS  AND  RECOMMENDATIONS 

1.  Unique  force  and  moment  data  are  presented  for  low  Froude 
number  operation  of  near-surface  hydrofoils. 

2.  At  shallow  submergences  (f/c  ~ 0.25)  there  are  remarkable 

changes  in  the  character  of  forces  and  moments  that  occur  in  the 

range  of  F = 0.5  to  0.75.  These  are  associated  with  a spanwide 
c 

hydraulic  Jump  or  breaking  wave  that  occurs  over  the  foil,  which 
in  turn  appears  to  be  caused  by  foil  thickness  as  well  as  angle  of 
attack.  These  changes  may  be  of  more  than  academic  interest  when 
consideration  is  given  to  schemes  for  partial  hydrofoil  support 
systems  of  large  ships. 

3.  The  measured  values  of  residual  drag  contribute  significantly 
to  the  total  drag  at  low  Froude  number  on  an  aspect  ratio  4 foil  at 
constant  lift,  as  determined  from  an  example  estimation  for  an 
18-foot-chord  hydrofoil  operating  at  f/c  = 1. 

4.  The  variation  of  hydrofoil  drag  at  low  Froude  numbers 
apparently  Includes  Important  contributions  due  to  the  presence  of 

a strut  and  pod.  These  should  be  better  understood,  both  experimentally 
and  theoretically. 
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5.  Additional  low  Froude  number  experiments  should  be  conducted 
with  forward-leading  sting (s)  support  in  order  to  better  Isolate  the 
lifting  planform  from  interference  drag  effects  caused  by  the  strut 
and  pod . 

6.  A cambered  foil  shape  should  be  used  in  order  to  achieve 
maximum  lift  coefficients  on  the  order  of  0.8  or  higher.  The  planform 
should  have  an  aspect  ratio  of  5 or  6,  and  some  consideration  should 
be  given  to  other  aspect  ratios  as  well. 

7.  Negative  angles  of  attack  should  be  included  in  any  future 
test  program. 

8.  Analytical  methods  and  computer  programs  to  implement  them 
should  be  developed  for  the  prediction  of  wavemaking  drag  on  hydrofoil- 
pod-strut  systems  at  arbitrary  Froude  number.  These  would  be  useful 
for  the  extraction  of  foil-only  data  from  experiments  and  for 
possible  optimization  of  very  large  hydrofoil  support  systems. 
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ASPECT  RATIO  A = 4 
SUBMERCENCE  RATIO  f/c  = 0.59 


Figure  1 - Predicted  Behavior  for  Lift  and  Invlscld  Drag 
Coefficients  Determined  by  Nlshlyama^ 
for  Rectangular  Planform  Hydrofoil 
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STRUT 


and  Drag  on  a Fully  Submerged  Hydro 


2:1  ELLIPTIC 


Geometry  of  the  Shroud  and  Strut 
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Figure  4 - Sketch  of  Experimental  Setup  Mounted  to  Carriage  II 
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Figure  9 - Measured  Drag  and  Lift  Coefficients  Versus  Chord 
Froude  Number  for  Submergence  Ratio  f/c  = 1.5 
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Figure  11  - Measured  Drag  and  Lift  Coefficients  Versus  Chord 
Froude  Number  for  Submergence  Ratio  f/c  = 3.5 
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Figure  15  - Measured  Drag  and  Lift  Coefficients  Versus  Submergence 

Ratio  for  Chord  Froude  Number  F =2.0 
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Measured  Drag  and  Lift  Coefficients  Versus  Angle 
of  Attack  for  Submergence  Ratio  f/c  * 0.25 
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Figure  21  - Measured  Drag  and  Lift  Coefficients  Versus  Angle 
of  Attack  for  Submergence  Ratio  f/c  -2.0 
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Figure  22  - Measured  Pitching  Moment  Coefficients  and  Center 
of  Pressure  Ratios  Versus  Chord  Froude  Number  for 
Submergence  Ratio  f/c  " 0.25 
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Figure  24  - Measured  Pitching  Moment  Coefficients  and  Center 
of  Pressure  Ratios  Versus  Chord  Froude  Number  for 
Submergence  Ratio  f/c  » 0.75 
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Figure  26  - Measured  Pitching  Moment  Coefficients  and  Center 
of  Pressure  Ratios  Versus  Chord  Froude  Number  for 
Submergence  Ratio  f/c  “ 1.5 
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Figure  27  - Measured  Pitching  Moment  Coefficients  and  Center 
of  Pressure  Ratios  Versus  Chord  Froude  Number  for 
Submergence  Ratio  f/c  » 2.0 
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Figure  28  - Measured  Pitching  Moment  Coefficients  and  Center 
of  Pressure  Ratios  Versus  Chord  Froude  Number  for 
Submergence  Ratio  f/c  - 3.5 
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Figure  29  - Lift-Slope  Coefficient  Versus  Chord  Froude 
Number  for  Various  Submergence  Ratios 
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Figure  30  - Relative  Lift-Slope  Coefficients  Versus 
Submergence  Ratios  for  Several  Chord 
Froude  Numbers 
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Figure  31  - Comparison  of  Lift  Coefficient  Versus  Angle  of  Attack 
for  Constant  Submergence  Ratio  and  Chord  Froude  Number 
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Figure  34  - Comparison  of  Lift  Coefficient  Versus  Angle  of  Attack 
for  Constant  Submergence  Ratio  and  Chord  Froude  Number 
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Figure  42  - Comparison  of  Theoretical  Inviscid  Drag  Coefficient 
Ratio  with  Experimental  Residual  Drag  Coefficient 
Ratio  for  a Rectangular  Hydrofoil 


Figure  44  - Residual  Drag  Coefficient  Versus  Lift  Coefficient 
Squared  for  Submergence  Ratio  f/c  =0.5 


Residual  Drag  Coefficient  Versus  Lift  Coefficient 
Squared  for  Submergence  Ratio  f/c  • 0.75 


Figure  46  - Residual  Drag  Coefficient  Versus  Lift  Coefficient 
Squared  for  Submergence  Ratio  f/c  =>  1.0 
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APPENDIX  B 


ZERO  AND  INFINITE  FROUDE  NUMBER  LIMITS  OF  HYDROFOIL  DRAG  DUE-TO-LIFT 

The  expressions  for  the  drag  coefficient  due-to-llft  in  the 
limits  of  infinite  and  zero  Froude  number  given  in  Equations  (8)  and 
(9),  respectively,  contain  the  biplane  function  o(X).  This  function 
depends  solely  on  the  parameter  X which  for  rectangular  planforms  is 


A formula  for  o(X)  valid  for  elliptic  foil  loading  can  be  deduced 

2 

from  results  discussed  by  Wu 


o(X)  * 1 - ^ / 1 + X2  [K(k)  - E(k)] 

— 1/2 

where  k = (1  + X^) 


K(k)  and  E(k)  * complete  elliptic  Integral  of  first  and 
second  kind,  respectively 

Approximating  formulas  for  K(k)  and  E(k)  are  given  by  Dwight^^  and 

2 

have  been  used  to  obtain  the  table  values  of  /C,  at  F -►0  and 

D^  L c 

F -»■<«>. 
c 

ZERO  AND  INFINITE  FROUDE  NUMBER  LIMITS  OF  C_  /C  ^ 

U ^ L> 


Zero  F 

Infinite  F 

c 

c 

f /c 

X 

o(X) 

1 + 6-0 
ttA 

1+6  + 0 
ttA 

0,25 

0.125 

0.6 

0.0342 

0.1297 

0.5 

0.25 

0.42 

0.0485 

0.1154 

0.75 

0.375 

0.31 

0.0573 

0.1066 

1.0 

0.5 

0.232 

0.0635 

0.1004 

1.5 

0.75 

0.14 

0.0708 

0.0931 

2.0 

1.0 

0.091 

0.0747 

0.0892 

3.5 

1.75 

0.0393 

0.0788 

0.0851 
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